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ABSTRACT 

In this paper, we investigate the relevance of using the 12 CO line emissions as in- 
dicators of star formation rates (SFR). For the first time, we present this study for 
a relatively large number of 12 CO transitions (12) as well as over a large interval 
in redshift (from z^O to z^6). For the nearby sources (D^IO Mpc), we have used 
homogeneous sample of 12 CO data provided by Bayet et al. (2004, 2006), mixing 
observational and modelled line intensities. For higher-z sources (z 1), we have col- 
lected 12 CO observations from various papers and have completed the data set of line 
intensities with model predictions which we also present in this paper. Finally, for 
increasing the statistics, we have included recent 12 CO(1-0) and 12 CO(3-2) observa- 
tions of intcrmcdiate-z sources. Linear regressions have been calculated for identifying 
the tightest SFR- 12 CO line luminosity relationships. We show that the total 12 CO, 
the 12 CO(5-4), the 12 CO(6-5) and the 12 CO(7-6) luminosities are the best indicators 
of SFR (as measured by the far-infrared luminosity). Comparisons with theoretical 
approaches from Krumholz and Thompson (2007) and Narayanan et al. (2008) are 
also performed in this paper. Although in general agreement, the predictions made 
by these authors and the observational results we present here show small and inter- 
esting discrepancies. In particular, the slope of the linear regressions, for J upP er ^ 4 
12 CO lines are not similar between theoretical studies and observations. On one hand, 
a larger high-J 12 CO data set of observations might help to better agree with mod- 
els, increasing the statistics. On the other hand, theoretical studies extended to high 
redshift sources might also reduce such discrepancies. 

Key words: Galaxies: starburst-nuclei-ISM - Submillimeter - ISM: molecules - 
infrared: galaxies 



1 INTRODUCTION 

Fifty years ago, it has been showed that the star formation 
rate (hereafter SFR) is intimately linked with the reservoir 
of the gas from which stars are forming (Schmidt 1959). The 
Kennicutt-Schmidt power law parameterized on local galax- 
ies by Kennicutt (1998a, b); Kennicutt et al. (2007) has led 
to an SFR index of iV=1.4±0.15. Since late 90s, researchers 
have converted the Kennicutt-Schmidt law into a more in- 
teresting relationship connecting the SFR (traced by the 
infrared luminosity - hereafter R r ) to the mass of molecular 
gas, investigating various molecular tracers. Firstly, Sanders 
et al. (1991); Sanders & Mirabel (1996) have showed that 
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the SFR is roughly proportional to the 12 CO(J=1-0) line 
(hereafter CO(1-0)) luminosity (slope of 1.4-1.6, consistent 
with the Kennicutt-Schmidt law index). This SFR-CO(l- 
0) relationship has been broadly interpreted as an increase 
of star formation efficiency as a function of molecular gas 
mass (and density). More recent analysis focussing on trac- 
ers of dense molecular gas such as the HCN(l-O) line (Gao 
& Solomon 2004a,b; Wu et al. 2005), the HCN(3-2) line 
(Bussmann et al. 2008) or the CO(3-2) line (Yao et al. 2003; 
Narayanan et al. 2005), have shown that these transitions 
are likely to be better indicators of star formation rate than 
the total H2 content (traced by the CO (1-0) luminosity). 
Indeed, the CO(1-0) line can be excited at rather low den- 
sities (~ 10 2 — 10 3 cm~ 3 ) and low temperature (~ 5K above 
ground) whereas the higher-J CO transitions and the HCN 
lines trace denser and warmer gas, more closely connected 
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with the stars in formation. However, recent modelling work 
(Krumholz & Thompson 2007; Narayanan et al. 2008) show 
that this situation is more complex than it appears, involv- 
ing in particular the values of the molecular line critical den- 
sities as compared with the mean gas density of the observed 
regions. 

So far, none of observational or modelling work in- 
vestigate the relationship between the SFR and molecular 
gas tracers both for several transitions (>2) of the same 
molecule (here CO) and over several factors of redshift (here 
z^sO-6). In this paper, the influence of the high- J CO lines 
is especially studied. This is motivated by the fact that, 
in addition to obtain much more accurate estimation of the 
CO line luminosities when adding the high-J transitions (see 
Bayet et al. 2004, 2006), more subtle effect may be revealed 
by these lines, when considering their respective critical den- 
sities. 

The knowledge of the global star formation history, and 
for individual galaxies, of their actual and past star forma- 
tion rate, is a key item for understanding galaxy evolution, 
and for comparing with state-of-the-art models. Statistical 
relations linking star formation properties with other galaxy 
global characteristics established by observing nearby galax- 
ies, have been already used to understand the processes rul- 
ing the large scale star forming activities (e.g. Malhotra et al. 
2001; Boselli et al. 2002). Such trends can afterwards be used 
in the large scale cosmological models, which lack the spatial 
and temporal resolution to describe the local details of star 
formation in individual galaxies. 

In this paper, we have thus compared the CO lines emis- 
sions (from J=l-0 to J=12-ll transitions) with the infrared 
emission, for the nearby galaxies we have surveyed using 
the Caltech Submillimeter Observatory (CSO) (see Bayet 
et al. 2004, 2006). We extended this comparison to higher 
redshifted sources (z = 1.4-6.4) for which CO line emissions 
have been previously observed (Cox et al. 2002; Bertoldi 
et al. 2003; Pety et al. 2004; Walter et al. 2004; Greve et al. 
2005; Solomon & Vanden Bout 2005; Tacconi et al. 2006; 
Weifi et al. 2007). For these distant objects, we derived the 
total and the individual CO line emissions (from J=l-0 to 
J=12-ll) by the same approach as the one presented in 
Bayet et al. (2004, 2006). In order to strengthen the results 
in a more statistical point of view, we finally have included 
additional literature CO(1-0) and CO(3-2) data from Gao & 
Solomon (2004a, b); Yao et al. (2003) and Narayanan et al. 
(2005), respectively. 

The paper is divided straightforwardly as followed: Sect. 
2 presents the data sample while Sect. 3 lists the results we 
have obtained, focussing especially on the comparison with 
the model predictions from Krumholz & Thompson (2007) 
and Narayanan et al. (2008). Finally, we conclude in Sect. 
4. 



2 SAMPLE SELECTION 

It has been crucial to determine, or find in the literature, the 
total CO line emission as well as the individual line emis- 
sions of the CO transitions from J=l-0 to J=12-ll. The 
total infrared luminosity (from 8 to 1000 /jm) for both the 
nearby and the high-z sources has also been crucial to deter- 
mine. Here, we have converted all the gathered data, for the 



Table 1. List of the detected CO lines for nearby galaxies in- 
cluded in this study. These CO lines detections arc from Bayet 
et al. (2004, 2006). Similar information for distant sources is found 
in Fig. 1 (observations represented by black bullets). When the 
line is not detected, we have used its corresponding predicted 
emission from modelling work (see Subsects. 2.1 and 2.2). For the 
Antennae sources (NGC 4038 and Overlap), the 12 CO(4-3) and 
12 CO(7-6) lines have been observed by Bayet et al. (2006) but 
they suffer from large uncertainties. Therefore, we used for these 
two objects and two transitions, rather the corresponding best 
model predictions. 



Name 


Distance 


detected 


12 CO lines 




Nearby sources 


(Mpc) 










IC 10 


1.0 1 


1-0, 2-1, 


3-2, 4-3, 


6-5, 


7-6 


NGC 253 


2.5 2 


1-0, 2-1, 


3-2, 4-3, 


6-5, 


7-6 


IC 342 


1.8 3 


1-0, 2-1, 


3-2, 4-3, 


6-5, 


7-6 


He 2-10 


9.0 4 


1-0, 2-1, 


3-2, 4-3, 


6-5, 


7-6 


NGC 4038 


13. 8 5 


1-0, 2-1, 


3-2, 4-3, 


6-5, 


7-6 


Overlap" 


13. 8 5 


1-0, 2-1, 


3-2, 4-3, 


6-5, 


7-6 


M 83 


3.5 6 


1-0, 2-1, 


3-2, 4-3, 


6-5 




NGC 6946 


5.5 7 


1-0, 2-1, 


3-2, 4-3, 


6-5 





°: Overlap corresponds to a shifted position from the NGC 4039 
nucleus in the Antennae Galaxy which shows a high gas density 
(see Bayet et al. 2006). References: 1 : Massey & Armandroff 
(1995); 2 : Adopted value from Mauersberger et al. (1996); 3 : 
McCall (1989); 4 : Kobulnicky & Johnson (1999); 5 : Savianc 
et al. (2004); 6 : Thim et al. (2003) and 7 : Tully (1988). 

first time, into a coherent and homogeneous sample we have 
corrected for various effects, as described below. In Table 1 
and Fig. 1, the observed CO data for nearby and distant 
sources are presented, respectively. Model predictions go up 
to the J=15-14 transition of CO, consistently to the work of 
Bayet et al. (2004, 2006). However, we have restricted our 
study to only the first 12 CO transitions, higher-J CO lines 
having very weak intensity (see Bayet et al. 2004, 2006 and 
Fig. 1). In the rest of the paper, the total CO will thus refer 
to the sum of only the first 12 CO transitions. 



2.1 Nearby sources 

The CO data we are using in this study for the nearby 
sources are from the observational and modelling work of 
Bayet et al. (2004, 2006). They already provided a consis- 
tent line intensity sample for 8 nearby (<10 Mpc) galaxies: 
NGC 253, IC 10, IC 342, NGC 6946, M 83, Henize 2-10 
and the Antennae (NGC 4038 and Overlap) from CO(1-0) 
to CO(15-14). We converted the first twelve CO integrated 
line intensities into luminosities (in Kkms _1 pc 2 ) using the 
formulae in Solomon & Vanden Bout (2005): 

L'co = 3.25 x 10 7 x Sco^v x v~ h \ x (1 + z)" 3 x D\ (1) 

where z is the redshift 1 , v \, s is the observed frequency 2 
(in GHz) of the CO line, ScoAu is the CO integrated line 
intensities (in Jykms -1 ) and Dl is the luminosity distance 



1 For nearby sources we used the NED database values of z. 

2 The observed frequency is equal to the rest frequency u rest 
divided by (1+z). 
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(in Mpc) . The total CO luminosity has been obtained by 
summing the individual CO line luminosities converted pre- 
viously from Eq. 1 into solar luminosity units (Lq) to avoid 
any frequency bias in the SFR-total CO relationship. Table 
2 summaries the results obtained. 

The infrared data used for the nearby sources and, more 
generally in the whole paper, are from Sanders et al. (2003). 
For consistency with the additional literature data described 
in Subsect. 2.3, we chose to use the total infrared luminosity 
from 8 fim to 1000 /im (LiR-.s-ioooum), thereby including 
the continuum emission shorter than 60 /j,m affected by the 
contribution from very small grains. This total infrared lu- 
minosity is known to be contaminated by a possible AGN 
contribution (see also Gracia-Carpio et al. 2008). In partic- 
ular, the AGN yield can be fairly large in the MIR range 
(e.g Rowan-Robinson & Crawford 1989) in active galaxies 
such as LIRGs and ULIRGs. However, for the sample of the 
nearby sources we are studying here, this AGN contamina- 
tion is considered as negligible since none of the 8 sources 
(IC 10, NGC 253, IC 342, Henize 2-10, NGC 4038, Overlap, 
M 83 and NGC 6946) is known for hosting an AGN (see 
e.g. Leroy et al. 2006; Martin et al. 2006; Usero et al. 2006; 
Wilson et al. 2000; Kramer et al. 2005; Israel & Baas 2001, 
respectively) . 

A more important aspect for nearby sources is that 
the CO measurements (both observations and model predic- 
tions) referred to an aperture of 22", as described by Bayet 
et al. (2004, 2006), while the infrared data from Sanders 
ct al. (2003) showed a higher aperture (80"). Thus, in Fig. 2 
where we present various SFR-CO luminosity relationships 
(see Sect. 3), we compare, for nearby sources, two emit- 
ting regions which are not spatially identical. Indeed, the 
region which emits the infrared luminosity is larger than 
the region where CO is detected. It is the same case for the 
intermediate-z data (Subsect. 2.3) but this problem does 
not appear in the case of higher-z sources since they are 
seen as point-like sources in both the sub-mm/mm and in- 
frared wavelengths. To correct this effect on nearby and 
intermediate-z sources, we have estimated a factor between 
the CO and the infrared luminosities resolution via the dust 
emission traced by SCUBA maps at 850 (im. We have de- 
rived from these maps the luminosity (removing the back- 
ground contribution) of the 850 /im emission at a spatial 
resolution of 22" and 80", for the sources available in the 
SCUBA archive. We have then calculated the ratio between 
the emissions observed at 22" and 80" that we have applied 
to the infrared data. We have obtained ratios varying from 
2.3 to 6.5, depending on the source. Rather than the 450 
/im SCUBA maps, we have chosen the 850/im SCUBA maps 
because they show a better signal-to-noise ratio. Neverthe- 
less, we checked that the factors obtained at 850/im were 
in agreement with the ones at 450/xm (difference obtained 
being less than 2%). After having applied such correction to 
the infrared data of the galaxies in common to this paper 
and the SCUBA archive, we have however obtained similar 
results to those listed in Table 3 and presented in Fig. 2 



3 For all the sources, we obtained the lumi- 
nosity distance using the web site calculator of 
http://www.astro.ucla.edu/~wright/CosmoCalc.html, within a 
cosmology of H D = 77 kms _1 Mpc _1 , Cl M =0.27 and fly =0.73. 



which show non-aperture corrected data. The slopes did not 
show changes greater than 5% in their values, depending on 
the SFR-CO line luminosity relationship studied. Due to the 
fact that all the data could not be corrected consistently for 
this effect since the galaxies presented here have not been all 
observed by SCUBA (e.g IC 342), we thus have decided to 
keep non-corrected infrared data in Fig. 2, increasing how- 
ever the error bars on the slope values by 5% in Table 3. 

Both the (non-aperture corrected) infrared and the CO 
luminosities used in the Fig. 2 are listed in Table 2. 



2.2 High-z sources 

So far, a complete observed CO SED does not exist for J=l-0 
to J=12-ll neither for nearby nor high-z sources. 

At high redshift, various CO transitions have been how- 
ever already detected (e.g. Cox et al. 2002; Bertoldi et al. 
2003; Pety et al. 2004; Walter et al. 2004; Greve et al. 
2005; Solomon & Vanden Bout 2005; Tacconi et al. 2006; 
Weifi et al. 2007). Most of these data are summarized in 
Solomon & Vanden Bout (2005). In our study, we restricted 
the number of studied sources presented in Solomon & Van- 
den Bout (2005) to ten objects : 4C60.07, APM 08279 4 , 
Cloverleaf QSO, SMM J14011, VCV J1409, PSSJ2322, TN 
J0924, SDSS J1148, IRAS F10214 and HR 10, keeping only 
those which show at least two CO line detections. Indeed, 
to correctly constrain the models which estimate the miss- 
ing line emissions (LVG models, see Appendix A) and derive 
relevant predicted CO line intensities, it is crucial to have 
as many CO line detections as possible. This is why we have 
rejected from our study high-z sources presented in Solomon 
& Vanden Bout (2005) with only one CO line observed. 

To estimate the total CO luminosity as well as the in- 
dividual missing molecular CO line luminosities in these 
sources, we have used a single-component Large Velocity 
Gradient (LVG) model as done in Bayet et al. (2004, 2006). 
The obtained line intensity predictions (see Fig. 1) corre- 
spond to a beam size of 22", consistently with other data 
sets. The goal of the paper is to investigate the SFR- 
molecular CO line luminosity relationships. Thus, we will 
not discuss further the physical properties of the molecular 
gas we have obtained using the LVG model. Nonetheless, we 
present them in Appendix A. 

In Fig. 1, we have superimposed on the CO observa- 
tions (black bullets with error bars), the predicted emis- 
sions of both the observed and the missing CO lines (grey 
filled triangles). Similarly to the nearby sources case, the 
total CO luminosity has been obtained by summing the in- 
dividual CO line luminosities converted previously from Eq. 
1 into solar luminosity units (L©) to avoid any frequency 
bias in the SFR-total CO relationship. More precisely, we 
have used the observed velocity-integrated CO line fluxes 
(Sco Av in Jykms -1 ) when detected and the predicted val- 
ues when not. The redshift values used in Eq. 1 are those 
presented in Solomon & Vanden Bout (2005). 

We did not need to correct the CO line luminosities 
for any beam dilution effect since the high-z galaxies are 



4 For this source, we rather used more recent CO data from WeiB 
ct al. (2007) 
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point-like sources within telescope beams and show unre- 
solved emissions whatever these wavelengths. However, some 
of these distant objects are lensed. Thus, we have applied 
the factor of lens magnification listed in Solomon & Vanden 
Bout (2005) to all the CO line luminosities. The underlying 
assumption of such a correction is that the gravitational lens 
magnifies similarly the emission from compact and warm re- 
gions usually traced by high-J CO lines such as the CO(7-6) 
line and the regions more extended normally traced by typ- 
ical low-J lines such as the CO(l-O) transition. In principle, 
this is a relevant assumption since the properties of the lens 
does not depend on the magnified source but only on the 
properties of the galaxies separated the observer from the 
studied source. 

For the infrared luminosity, we used the integrated val- 
ues ("LF/fl(int.)") listed in Table 1 of Solomon & Vanden 
Bout (2005), corresponding to infrared emission corrected 
for lens magnification. After having checked several refer- 
ences in this table, it appears that the "LF/jj(int.)" values 
are similar to the L/^s-ioooum presented in Sanders et al. 
(2003) (see the example of HR 10 for which Dey et al. 1999 
derived the infrared luminosity by fitting the dust SED on 
a rest wavelength range from 10 /im to 2 cm.). 

For these distant sources, we also present in Table 2, 
both the infrared and the CO luminosities we have used in 
Fig. 2. 

2.3 Intermediate-z data from the literature 

We have included in our study CO data from observations 
at intermediate redshift (within similar beam sizes) to bet- 
ter constrain the SFR-molecular CO line luminosity rela- 
tionships, and especially increase the statistics. These data 
are either from Gao & Solomon (2004a,b) (CO(1-0) tran- 
sition observed with a 22" beam size), or from Yao et al. 
(2003) and Narayanan et al. (2005) (CO(3-2) line observed 
with a 15" and 22" beam sizes, respectively). We have not 
included the CO(2-l) data reported in Rigopoulou et al. 
(1996) (beam size of 23") because they appear inconsistent 
with those presented in Yao et al. (2003); Gao & Solomon 
(2004a, b); Narayanan et al. (2005). We have included in our 
SFR-CO relationships the observed and modelling results 
dedicated to Markarian 231 (Papadopoulos et al. 2007), pro- 
viding CO line luminosities estimations up to the CO(10-9). 
For the infrared luminosity values, we have used the values 
reported in Sanders et al. (2003) and not those listed in each 
corresponding papers. 

The additional sources included here are seen to be 
mostly LIRGs and ULIRGs (at a distance > 10 Mpc). Ra- 
dio continuum maps (Condon et al. 1990), and the HI and 
SCUBA 850 ^m maps (Thompson et al. 2002), show that 
much of the emission of the objects is extended with re- 
spect to the beam size of 15". However, in recent high- 
resolution (2"-3") millimeter wave observations of seven 
LIRGs/ULIRGs, Bryant & Scoville (1999) find that nearly 
all of the detected CO(J=1-0) emission is concentrated 
within the central 1.6kpc in six of seven objects. Because it 
is unlikely that there will be significant high-J CO emission 
where there is no CO(1-0) detected, we assume that all of the 
emitting gas in these additional sources is confined within 
the same region. Consequently, we have not converted the 
14" beam data into that for a 22" beam. 



As for nearby galaxies, the infrared luminosities of these 
sources are expressed for a larger aperture than the one used 
for the CO data. As previously mentioned (see Subsect. 2.1), 
this effect on the SFR-CO line luminosity relationship has 
been investigated through SCUBA data. Unfortunately, as 
for the case of nearby galaxies, the intermediate-z sources we 
are using here have indeed not been all observed by SCUBA, 
therefore no correction on their infrared data has been per- 
formed directly. We have however taken fully into account 
the effect of different spatial resolutions on the slope val- 
ues listed in Table 3 (increased error bars as explained in 
Subsect. 2.1). 

In the samples presented in Yao et al. (2003); Gao & 
Solomon (2004a,b); Narayanan et al. (2005), some sources 
such as ARP 220 are in common. For removing any calibra- 
tion effects on the SFR-CO relationships, we have followed 
the recommendation of Narayanan et al. (2005) and applied 
their scaling factor equal to 0.26 for the CO(3-2) and to 0.45 
for the CO(1-0). They specified that this effect may indeed 
occur when various instruments (and thus various calibra- 
tion processes) are used. 

Note that the intermediate-z sources increase signifi- 
cantly the number of CO observations used in the correla- 
tions either in the SFR-CO(l-O) line luminosity or in the 
SFR-CO(3-2) line luminosity relationships (see triangles in 
Fig. 2 as well as the fourth column in Table 3). This strength- 
ens the corresponding results. 



3 RESULTS 

3.1 General arguments 

We present in Fig. 2, the SFR-CO line luminosity relation- 
ships we have obtained for the total CO line emission, the 
CO(1-0), CO(3-2) and CO(7-6) line emissions. In the same 
vein, we have obtained the SFR-CO line luminosity relation- 
ships for the other CO transitions: J=2-l and from J=4-3 
to J=12-ll. SFR-CO luminosity relationships for CO data 
higher than CO (7-6) are less relevant than the correlation in- 
volving lower-J CO lines. Indeed, most of the data in CO(8- 
7), CO(9-8), CO(10-9), CO(ll-lO) and CO(12-ll) are from 
modelled fit, similarly to the fits presented in Fig. 1. In such 
cases, correlations between the SFR and the CO luminosity 
are not constrained by many observational measures. Thus, 
we exclude these correlations from the following analysis. 

As one can see in Fig. 2, a proportionality exists be- 
tween the CO and infrared luminosities over a large range 
of redshift (galaxies from z~ up to z~ 6 sources). This 
characteristic has been previously seen (Sanders et al. 1991; 
Sanders & Mirabel 1996 and references therein) although 
this sample was more restricted in redshift than the one we 
present in this paper. It has been interpreted as an increas- 
ing star formation efficiency (SFE; SFR divided by M(H2)) 
as a function of molecular gas mass. Here, our data sample 
confirms this broad interpretation. 

To quantify better these relationships, we performed 
linear regressions (using the software xmgrace). The corre- 
sponding output parameters are listed in Table 3. Note that 
IC 10 has been excluded from these calculations. For the 
SFR-CO(3-2) relationship we have also excluded NGC 7817, 
as Narayanan et al. (2005) suggested. Indeed these two 
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Figure 1. Observed and predicted CO spectral energy distributions (Jykms 1 ) of the following high rcdshift sources (from top to bottom, 
and by increasing order of redshift): HR 10, IRAS F10214, The Cloverleaf QSO, SMM J14011, VCV J1409, 4C60.06, APM 08279, PSS 
J2322, TN J0924 and SDSS J1148 (see plots). Observations and their corresponding error bars are represented by black bullets (see 
references listed in Solomon & Vanden Bout 2005) while LVG predictions (sec Appendix A) are symbolized with grey triangles. To make 
the plots more easily readable, we have connected the CO LVG predictions (grey lines). One could notice that, except for the source 
HR 10, the position of the maximum of the SED is located at rather high-J (J U pper ^4). 
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Table 2. Infrared and CO lines luminosities of the sources in our sample (see squares symbols in Fig. 2). For distant galaxies, we also 
list in the last column the lens magnification factor we have found in Solomon & Vanden Bout (2005) and use in this study (see Fig. 
2). Intermediatc-z sources values are not shown in this table since already described in details in the references mentioned in the text. 
When an asterisk is seen, it means that the value has been estimated from observations only. 



l\Tn t > i o 


T, a 
I R 


u CO(tot) 


Ref. 


L'co(i— 0) 


Ref. 


T 7 

u CO(3— 2) 


Ref. 


L'cO(7-6) 


Ref. 






L© 


L© 




Kkms - 1 pc 2 




Kkms - 1 pc 2 




Kkms - 1 pc 2 






IC 10 


8.2xl0 9 


4.4xl0 2 


1 & 2 


5.6xl0 6 * 


1 & 2 


3.5 x10 s * 


1 & 2 


7.3xl0 3 * 


1 & 2 




NGC 253 


2.8xl0 1() 


3.3xl0 5 


2 & 3 


1.5xl0 9 * 


1 & 2 


1.4x10 s * 


1 & 2 


8.0x10 s * 


1 & 2 




IC 342 


1.5xl0 1() 


8.4xl0 3 


1 & 2 


1.6xl0 8 * 


1 & 2 


9.7xl0 6 * 


1 & 2 


1.2x10 s * 


1 & 2 




He 2-10 


6.2xl0 9 


3.4X10 4 


2 & 3 


2.1x10 s * 


1 & 2 


l.lxlO 7 * 


1 & 2 


1.0x10 s * 


1 & 2 




NGC 4038 


6.9xl0 1() 


5.5xl0 5 


1 & 2 


2.1xl0 9 * 


1 & 2 


1.9x10 s * 


1 & 2 


1.3xl0 7 * 


1 & 2 




Overlap 


6.9xl0 1() 


6.7x10 s 


1 & 2 


3.3xl0 9 * 


1 & 2 


2.9x10 s * 


1 & 2 


1.7xl0 7 * 


1 & 2 




M 83 


1.3xl0 1() 


1.0x10 s 


1 & 2 


3.1x10 s * 


1 & 2 


3.3xl0 7 * 


1 & 2 


2.7x10 s 


1 & 2 




NGC 6946 


1.5xl0 10 


1.0x10 s 


1 & 2 


1.4xl0 9 * 


1 & 2 


1.1x10 s * 


1 & 2 


1.3x10 s 


1 & 2 




Name 


T a 


T ,b 

^ CO(tot) 


Ref. 


L'co(i-o) 


Ref. 


L 'cO(3-2) 


Ref. 


L'cO(7-6) 


Ref. 


Lens 




L© 


L© 




Kkms - 1 pc 2 




Kkms - x pc 2 




Kkms - 1 pc 2 




magn. factor 


HR10 


6.5xl0 12 


4.2xl0 8 


2 


5.6xl0 10 * 


4 & 2 


2.5X10 10 


2 


7.9x10 s 


2 


1 


IRAS F 102 14 


3.6xl0 12 


6.6xl0 9 


2 


1.4xlO n 


2 


9.9xl0 10 * 


4 & 2 


3.2xl0 10 * 


4 & 2 


17 


Cloverleaf 


5.4xl0 12 


4.3xl0 1() 


2 


5.0xlO n 


2 


3.9X10 11 * 


4 & 2 


1.9X10 11 * 


4 & 2 


11 


SMM J14011 


2.4xl0 12 


6.1xl0 9 


2 


l.OxlO 11 


2 


8.3xl0 10 * 


4 & 2 


1.7xl0 10 * 


4 & 2 


8.3 


VCV J1409 


3.5xl0 13 


7.8xl0 9 


2 


9.6xl0 9 


2 


6.9xl0 10 


4 & 2 


2.3xl0 10 


4 & 2 


1 


4CO60.07 


1.3xl0 13 


4.8xl0 9 


2 


6.1xl0 10 * 


4 & 2 


4.9xl0 10 


2 


l.OxlO 10 


2 


1 


APM 08279 


2.9xl0 13 


l.lxlO 11 


2 


8.0xl0 10 * 


5 & 2 


1.3X10 11 


5 & 2 


l.OxlO 11 


5 & 2 


7 


PSS J2322 


9.3xl0 12 


9.6xl0 9 


2 


l.lxlO 11 * 


4 & 2 


1.8X10 11 


2 


8.3xl0 9 


2 


2.5 


TN J0924 


7.0xl0 12 


9.9xl0 9 


2 


7.3xl0 10 * 


4 & 2 


6.0xl0 10 


2 


1.7xl0 10 


2 


1 


SDSS J1148 


2.7xl0 13 


9.1xl0 9 


2 


3.5xl0 9 


2 


2.2xl0 10 


4 & 2 


1.5X10 10 


4 & 2 


1 



a : Lm is the 8 - 1000 /im total infrared luminosity deduced from cither Sanders ct al. (2003) (nearby sources) or Solomon & Vanden 
Bout (2005) (high-z galaxies). For the distant objects, Lir are corrected from the lens magnification effect (see text in Subsect. 2.2); b : 
L 'co(tot) ' s the total CO line luminosity obtained after the use of models for predicting the emission of the missing (non-observed) CO 
lines (see the Subsccts. 2.1 and 2.2 and, for more details, the Appendix A and the modelling work of Bayet ct al. 2004, 2006). 
References : 1: See Bayet et al. (2006); 2: This work; 3: Bayet ct al. (2004); 4: Solomon & Vanden Bout (2005); 5: Weifi et al. (2007) 



sources show very small \^co value with respect to the lo- 
cus of other objects. This is especially true for IC 10 over 
several CO transitions (see Fig. 2). 

Deviations from the SFR-CO trend can be expected for 
subsolar metallicity such as IC 10 (12 log(0/H)= 8.31 ± 
0.2 from Arimoto et al. 1996). As shown by Lequeux et al. 
(1994), a CO deficit (and a [CII] excess), relatively to the in- 
frared luminosity is expected for low metallicity systems. At 
low metallicity, molecular hydrogen can still be formed on 
grains, although the threshold between atomic and molec- 
ular hydrogen is slightly shifted compared to what is seen 
in our Galaxy (see the discussion on the FUSE results by 
Tumlinson ct al. 2002). In addition, CO is more easily de- 
stroyed by FUV photons than H2. In diffuse and translucent 
clouds, the decrease of CO photodissociation with increasing 
A„ is provided by the combined effect of self-shielding and 
dust extinction. At lower metallicities, these two effects get 
weaker, hence CO can not survive at small A„. In molecular 
clouds, CO survives in the core, but an extended envelope 
where carbon is ionized is also present. Therefore, the CO 
emission is decreased, (and the [CII] emission increased) rel- 
ative to normal metallicity conditions. 

Over the last years, in the literature, other tracers of 
molecular gas have been investigated as SFR indicators, such 
as the HCN (Gao & Solomon 2004a,b; Vanden Bout et al. 
2004; Carilli et al. 2005; Bussmann et al. 2008). The linear 
and tighter relationship shown between the HCN(l-O) line 
luminosity and the infrared luminosity, as compared to the 



relation SFR-CO (1-0), has been first interpreted by the fact 
that the HCN is a molecule tracing warmer and denser gas, 
more closely link with the stars in formation. Indeed, the 
CO(1-0) tends to emit from both dense core and more diffuse 
molecular filaments and cloud atmosphere (n cr ; t ~ 10 2 — 10 3 
cm -3 ) leading to a non-linear relationship, whereas the HCN 
is typically only thermalized in the dense cores of molecular 
clouds (n cr it ~ 10 s cm -3 ), allowing thus linearity. However, 
the fact that, in a similar sample of galaxies Yao et al. (2003); 
Narayanan et al. (2005) also found a linear SFR-CO (3-2) 
relationship provide evidence against the sole HCN-related 
chemistry explanation (The critical density of the CO (3-2) 
is n cri t ~ 10 4 cm -3 ) . 

More recent theoretical analysis (Krumholz & Thomp- 
son 2007; Narayanan et al. 2008) have shown that the rela- 
tionship between the SFR and molecular line luminosity in 
star-forming clouds (and thus, the value of the slope of the 
linear regressions) depends rather on how the critical density 
of the molecular transition compares to the mean density of 
the observed sources. Lines with critical densities smaller 
than the mean density in a observed region (e.g. CO(l- 
0)) probe the total molecular gas and the SFR-molccular 
line relationship slope is in agreement with the Kennicutt- 
Schmidt index (slope ~ 1.5). On the contrary, lines with 
critical densities larger than the mean density (e.g. HCN(1- 
0) and CO (3-2)) show luminosities rising faster than linearly 
with increasing mean gas density and the SFR-molecular 
line relationship is sub-linear (slope < 1.0). In this paper, 
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we confirm observationally these conclusions (see Table. 3) 
extending them, in addition, to a larger interval in redshift. 

3.2 More detailed analysis 

Fig. 2 and Table 3 show that the SFR -CO(5-4),-CO(6-5) 
and -CO(7-6) line luminosity relationships are the tightest 
correlations (highest values of the correlation coefficients). 
Thus, between all the CO transitions, for the data sample 
we present here, these transitions might be considered as 
the best star formation rate indicators in galaxies. The total 
CO luminosity is also a very good tracer of star formation 
rate in galaxies since it shows a tight correlation (correlation 
coefficient of 0.984). 

These tight SFR-high-J CO correlations contrast with 
the one seen for the CO(f-O) line. Indeed, even if we restrict 
the sample of CO(I-O) observations to only the 17 sources we 
used in the SFR-higher-J CO line luminosity relationships, 
considering like that the same statistics, the correlation co- 
efficient of the CO(I-O) line is still the lowest (see Table 
3). This has been previously mentioned by Riechers et al. 
(2006). In the same vein, it has been shown by Bussmann 
et al. (2008) that the higher-J HCN line such as the HCN(3- 
2) transition is more tightly correlated to the SFR traced by 
the infrared luminosity than the HCN(f-O) line. We could 
thus conclude that the CO(1-0) line is not a good indicator 
of star formation rates in galaxies. 

In a more detailed analysis of the plots seen in Fig. 2, 
it is interesting to note that, even if the linear regressions 
fit rather well in all cases the observations, they however 
become slightly not so relevant at low CO luminosity, and 
thus, especially for the SFR-high-J CO relationships. This 
tendency is even more obvious in the SFR-CO(3-2) plot. The 
galaxies IC 342 and Henize 2-fO seem to be mainly respon- 
sible of this divergence. In the SFR-CO(7-6) plot as well as 
in the correlations involving higher-J CO lines, it is more 
difficult to see this trend, probably because the statistics is 
reduced to f 7 sources. This leads to a less accurate linear re- 
gression than in the SFR-CO(3-2) case. In the SFR-CO(7-6) 
line luminosity relationship, one could also notice that the 
deficit of CO luminosity of the galaxy IC 342 is even more 
enhanced. The fact that IC 342 is a less active galaxy and 
that Henize 2- 10 is small in our nearby source sample may 
also play a role. Indeed, as compared to more active galaxies 
(such as NGC 253), they may contain less dense gas. This 
has been recently confirmed by Bayet et al. (2009) who de- 
tected the CS(7-6) line with much more difficulties in IC 342 
than NGC 253 (signal-to-noise ratio lower by a factor of 3 
between the two sources). Consequently more observations 
of high-J CO lines in quiet nearby galaxies are needed to 
explain better this characteristic. Another possible explana- 
tion may also be the degree of thermalization varying from 
one environment to another as explained in Narayanan et al. 
(2008). 

In the same way, it is noticeable that a small discrep- 
ancy between the linear regressions and the observations ap- 
pears for sources showing the highest infrared luminosities 
(see Fig. 2 in the upper right corners). However, contrarily to 
the previous case (nearby sources), this effect is not seen for 
the SFR-CO(7-6) plot nor for the other SFR-higher-J CO 
correlations. On the contrary, it appears clearer for lower-J 
CO lines. One very simple first assumption should be that 



this discrepancy may be caused by the fact that such ac- 
tive (high-z) galaxies is expected to have a higher percent 
of dense gas than in local galaxies. Indeed these objects are 
known to undergo more violent star-burst phases. The gas 
component traced by the CO(f-O) line may not be (rela- 
tively) as abundant and as bright as in the local Universe. 
This hypothesis may be confirmed by the localization of the 
CO SED maximum seen rather at high-J CO lines (see Fig. 
1) which are known to trace a warmer and denser gas (see 
Bayet et al. 2004, 2006). 



3.3 Comparison with model predictions 

It is essential to confront the results we have obtained to 
the model developed by Krumholz & Thompson (2007); 
Narayanan et al. (2008). More especially, Narayanan et al. 
(2008) provided some estimations of linear regression slopes 
for various SFR-CO line luminosity relationships. In Fig. 
3, we have reproduced the Figure 7 of Narayanan et al. 
(2008) adding the new observational constraints we have 
obtained (see the open white circles with error bars in Fig. 
3). Fig. 3 thus provides the variation of the slope of the 
SFR-CO luminosity relationship with respect to the transi- 
tion of CO investigated. The best predicted (from models) 
SFR-CO slopes are represented by a dashed black line while 
the horizontal grey lines symbolize its scatter. This scatter 
is computed by randomly drawing a sample of 19 galaxies 
(which is similar to the size of our nearby and high-z source 
sample) out of a set of 100 model galaxies 1000 times. Here, 
our measurements are consistent with other measurements 
(see the crosses in the Fig. 3) and with model predictions 
up to the SFR-CO(3-2) relationship. However, for higher- 
J CO line, some discrepancies are seen. As mentioned in 
Narayanan et al. (2008), these discrepancies may likely due 
to the fact that these models do not include any high-z 
sources. Most of these sources have an increasing fraction 
of AGN as compared with this of the local galaxies. These 
embedded AGN might thus significantly contaminate the in- 
frared luminosity of the sources (Tran et al. 2001; Kim et al. 
2002; Veilleux et al. 2002), as well as be responsible of an 
increase of the gas temperature, leading to a higher slope 
in the linear regression. In a simple view, the star-formation 
being warmer in such distant objects, the temperature effect 
is expected especially enhanced for the high-J CO line re- 
lationships. Temperature contamination, imperfectly taken 
into account in the models, is thus more expected in such 
CO lines than in low-J CO transitions. 



4 CONCLUSIONS 

In this paper we have investigated the relevance of using 
the CO line emissions as indicators of star formation rate. 
For the first time, we have both studied this question over a 
large number of CO transitions (12) and over a large sample 
of source going up to a redshift of z>6. We have shown that 
the SFR-total CO, -CO (5-4), -CO (6-5) and -CO (7-6) line lu- 
minosity relationships are the tightest correlations, making 
these lines the best indicators of the star formation rates. 
The results we have obtained also strongly confirmed the 
predictions from Krumholz & Thompson (2007); Narayanan 
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Figure 2. Examples of SFR-CO line luminosity relationships we have obtained, expressed in log-log scales. From left to right, we plotted 
the total integrated infrared luminosity L/i{ : g_iooOM"i ( m Lq) versus the total CO (in Lq), the CO(l-O), the CO(3-2) and the CO(7-6) 
luminosities (in Kkms - ^c 2 ) (see Sect. 2). The squares correspond to the sources we studied in details in this paper (see Subsects. 2.1 
and 2.2) while the triangles correspond to the literature data we have added from Yao ct al. (2003); Gao & Solomon (2004a, b); Narayanan 
et al. (2005) (see Subscct. 2.3). The black filled squares represent the nearby sources while the white opened squares symbolize the high-z 
sources. In all plots, IC 10 is separated from other nearby sources (see Subsect. 3.1) and represented by open white circles. In the plot 
representing the SFR-CO(3-2) luminosity relationship one can see a lambda symbol corresponding to the source NGC 7817 we have also 
excluded (see Subscct. 3.1). The linear regressions have been obtained using the xmgrace software. They are represented by black lines. 
For keeping the figures clear, we plotted the typical errors bars of the observations (included in the calculations of the linear regressions), 
only on the SFR-CO (7-6) luminosity correlation. These errors correspond to a ±lcr uncertainties both on the Ljjj and on the molecular 
CO line luminosities. 



Table 3. Results of the linear regressions (slope and correlation 
coefficients) for the relationships between SFR and the CO tran- 
sitions from J=l-0 to J=12-ll (see Fig. 2 for some examples of 
SFR-CO line correlations). The linear regressions have been ob- 
tained using the software xmgrace, including the error bars of 
the observations in the calculations. In these linear regressions, 
we excluded two sources from the calculations: IC 10 for all the 
regressions and NGC 7817 for the SFR-CO(3-2) relationship (see 
Sect. 3.1). We have isolated the SFR- total CO luminosity rela- 
tionship from others because we have used different units. 



Molecular line 


Slope 


Correlation 


Number 


Luminosities 




Coeff. 


of sources 


CO(l-0) a 


1.41±0.26 


0.82 


17 


coii-o) 6 


1.14±0.08 


0.82 


103 


CO(2-l) 


1.20±0.10 


0.95 


17 


CO(3-2) a 


1.00±0.07 


0.97 


17 


CO(3-2) 6 


0.99±0.06 


0.89 


81 


CO(4-3) 


0.94±0.05 


0.98 


17 


CO(5-4) 


0.90±0.04 


0.98 


17 


CO(6-5) 


0.86±0.04 


0.98 


17 


CO(7-6) 


0.80±0.04 


0.98 


17 


CO(8-7) 


0.74±0.05 


0.97 


17 


CO(9-8) 


0.67±0.06 


0.94 


17 


CO(10-9) 


0.61±0.07 


0.92 


17 


CO(ll-lO) 


0.57±0.07 


0.90 


17 


CO(12-ll) 


0.53±0.07 


0.89 


17 


CO(tot) 


0.62±0.04 


0.97 


16 c 



a : Without any additional literature data; b : With additional lit- 
erature data; c : Markarian 231 is not included in this linear re- 
gression calculation because we did not find any estimation of its 
CO(ll-lO) and CO(12-ll) lines luminosities needed for estimat- 
ing properly its total CO luminosity. 



2.0 1 1 1 1 1 1 r 



1.5 - 




0.0 1 i i i i i i i 

CO1-0 CD2-1 C03-2 C04-3 C05-4 C06-5 C07-6 
Transition 

Figure 3. Predicted slopes in log(L/^)-log(L' C 7o) space as a 
function of J transitions of CO, as seen in the study of Narayanan 
ct al. (2008) (Figure 7). The grey horizontal lines are the model 
predictions and the crosses, some observations Narayanan ct al. 
(2008) used. Observational constraints on the slope measured 
from CO(1-0) to CO(7-6) presented in this paper are symbol- 
ized by open white circles (with error bars). From the CO(4-3) 
line, the best-fit slopes (dashed line) given by the predictions from 
Narayanan et al. (2008) become not consistent with our observa- 
tional results. 



et al. (2008) who showed that the SFR-CO lines luminos- 
ity relationships are above all regulated by the Kennicutt- 
Schmidt law, which sets the way in which observed transi- 
tions trace molecular gas. In other words, the SFR-CO line 
luminosity relationships depends on how the critical den- 
sity of the molecule compares to the mean density of ob- 
served source. We confirm in this paper that, for the CO 
lines with J upP er >3, the SFR-CO line luminosity relation- 
ship are indeed sub-linear. However some discrepancies be- 
tween model predictions and observations appear for higher- 
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J CO lines (from CO(4-3)) relationships. These differences 
may be likely due to the fact that, in the model the high-z 
sources are not included (possible temperature contamina- 
tion from AGN heating processes). Anyway, more observa- 
tions (ALMA, Herschel) of CO in both nearby and high-z 
sources, especially the high frequency CO transitions could 
also be very helpful for explaining these discrepancies. 
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APPENDIX A: SINGLE COMPONENT LVG 
MODELLING RESULTS FOR THE HIGH-Z 
SOURCES 

As mentioned in Subsect. 2.2, a complete observed CO SED 
does not exist for J=l-0 to J=12-ll in high-z sources. We 
thus have collected CO line emissions from Solomon & Van- 
den Bout (2005) and performed a Large Velocity Gradient 
(LVG) analysis to derive the missing CO integrated line in- 
tensities. The results of this modelling work is shown in Fig. 
1 and in the following Table Al. 

We have run LVG models (Goldreich & Kwan 1974; de 
Jong et al. 1975), already well described in various papers. 
The version we have used in this study is the one presented 
in Bayet et al. (2004, 2006, 2008). This LVG model is ba- 
sically running with three free parameters: the gas density 
n(H_2), the kinetic temperature (T*,) and the CO column 
density divided by the line width (N( 12 CO)/Av). We have 
investigated the following range of LVG input parameters : 
5 K <T K < 255 K, 1 x 10 12 cm-7kms _1 <N( 12 CO)/Ai; < 
1 x 10 19 cm" 2 and 1 x lOW" 3 <n(H 2 )< 1 x 10 7 cm~ 3 . 

For determining the missing CO line intensities, we have 
constrained the predicted CO line intensity ratios with the 
observed values (computed from Solomon & Vanden Bout 
2005) via a reduced % 2 method as performed in Bayet et al. 
(2004, 2006, 2008). 

We have used a single LVG component for reproduc- 
ing the entire set of CO data in high-z sources, being fully 
conscious that this modelling approach is not ideal for re- 
producing the widespread view that the ISM contains dense, 
star-forming molecular cloud cores and more diffuse gas in 
cloud envelopes, even in high-z sources. Motivated also by, 
on average, the small number of CO detections, we have 
thus considered this approach as reasonable for the purpose 
of this paper. 

We recommend to use the predicted physical properties 
derived and listed in Table Al as only indicative values. 
We indeed remind the reader that we were not trying to 
model each source individually, deriving the accurate set 
of physical properties for the molecular gas, but that we 
rather aim to obtain satisfactory estimations of the CO line 
intensities in these sources. Note that, for APM 08279 (Weifi 



Table Al. Results of the single component LVG model analysis. 
We presented in this table the physical properties of the best 
LVG model (having the lowest-x 2 value when compared to the 
observations) . 



Source 


T K 


n(H 2 ) 


N(CO)/Au 




in K 


( cm- 3 ) 


(cm _2 /kms _1 ) 


HR10 


255 


l.OxlO 1 


8.0xl0 18 


IRAS F10214 


15 


8.0xl0 5 


3.0xl0 18 


Cloverleaf 


20 


4.0xl0 4 


l.OxlO 19 


SMM J14011 


205 


2.7X10 1 


7.9xl0 18 


VCV J1409 


185 


1.7xl0 4 


4.5X10 15 


4CO60.07 


180 


1.7X10 2 


1.5xl0 18 


APM 08279 


255 


9.0xl0 3 


8.5xl0 17 


PSS J2322 


30 


5.0xl0 4 


2.2xl0 16 


TN J0924 


205 


1.4xl0 2 


2.5xl0 18 


SDSS J1148 


100 


7.9xl0 4 


4.5xl0 15 



et al. 2007), the use of two LVG components has been made 
but, such as the sum of both CO components agrees with the 
observed CO SED. Thus, the choice of a one or a two LVG 
components model does not affect significantly the SFR-total 
CO luminosity relationships we present in this paper. 

This paper has been typeset from a TfrX/ DTjtX file prepared 
by the author. 



